Clustering Fossils from the Early Universe 
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Many inflationary theories introduce new scalar, vector, or tensor degrees ol freedom that may 
then affect the generation of primordial density perturbations. Here we show how to search a galaxy 
(or 21-cm) survey for the imprint of primordial scalar, vector, and tensor fields. These new fields 
induce local departures to an otherwise statistically isotropic two-point correlation function, or 
cquivalently, nontrivial four-point correlation functions (or trispectra, in Fourier space), that can be 
decomposed into scalar, vector, and tensor components. We write down the optimal estimators for 
these various components and show how the sensitivity to these modes depends on the galaxy-survey 
parameters. New probes of parity-violating early-Universe physics are also presented. 

PACS numbers: 98.80.-k 



INTRODUCTION 

Galaxy clustering has proven to be invaluable in as- 
sembling our current picture of a Universe with a nearly 
scale-invariant spectrum of the primordial curvature per- 
turbations [lj. The principle tool in clustering stud- 
ies has been the two-point correlation function, or in 
Fourier space the power spectrum, determined under the 
assumption of statistical isotropy (SI). With the advent 
of new generations of galaxy surveys, as well as longer- 
term prospects for measuring the primordial mass distri- 
bution with 21-cm surveys of the epoch of reionization 
[2] and/or dark ages [5], it is worthwhile to think about 
what can be further done with these measurements. 

Many inflationary models introduce new fields that 
may couple to the inflaton responsible for generating cur- 
vature perturbations. The effects of these fields may then 
appear as local departures from SI, or as non-Gaussianity, 
in the curvature perturbation. For example, models with 
an additional scalar field introduce a nontrivial four-point 
correlation function (or trispectrum, in Fourier space) [4] , 
which we below will describe as local departures from sta- 
tistical isotropy; apart from this correlation, the scalar 
field may leave no visible trace. There may also be vec- 
tor (spin-1) fields U M [5] — or vector spacetime-metric per- 
turbations brought to life in alternative-gravity theories 
[B] — that, if coupled to the inflaton p (e.g., through a 
term (d^p>)(d u p)d tl V v ) may leave an imprint on the pri- 
mordial mass distribution without leaving any other ob- 
servable trace. Similar correlations with a tensor (i.e., 
spin-2) field T* v (e.g., {d^){d„ip)T^) can be envi- 
sioned. Even in the absence of new fields, there are ten- 
sor metric perturbations (gravitational waves) that may 
be correlated with the primordial curvature perturbation 
H]. Tensor distortions to the two-point correlation 
function ( "metric shear" ) may also be introduced at late 
times [9j [10] , and late-time nonlinear effects may induce 
scalar-like distortions to the two-point function [IT] . 

Here we describe how the fossils of primordial tensor, 



vector, and scalar fields are imprinted on the mass distri- 
bution in the Universe today. We express these relics in 
terms of two-point correlations that depart locally from 
SI or off-diagonal correlations of the density-field Fourier 
components. This formalism allows the correlations to be 
decomposed geometrically into scalar, vector, and tensor 
components. We write down the optimal estimators for 
these scalar, vector, and tensor correlations and quantify 
the amplitudes that can be detected if these perturba- 
tions have (as may be expected in inflationary models) 
nearly scale-invariant spectra. 

We begin with the null hypothesis that primordial den- 
sity perturbations are statistically isotropic and Gaus- 
sian. This implies that the Fourier modes <5(k) of the 
density perturbation <5(x) (at some fixed time) have co- 
variances, (<5(k)(5(k')) = V8P_ k ,P(k), where the Kro- 
necker (Dirac) delta on the right-hand side is zero unless 
k = — k', P{k) is the matter power spectrum, and V is 
the volume of the survey. In other words, the different 
Fourier modes of the density field are uncorrelated under 
the null hypothesis. 

Coupling of the inflaton to some other field produces 
non-Gaussianity in the mass distribution that appears as 
off-diagonal (i.e., ki ^ — k 2 ) correlations of the density- 
field Fourier components in the presence of a given real- 
ization of the new field. Global SI requires that a given 
Fourier mode h p (K.) of wavevector K and polarization p 
(about which we will say more below) of the new field 
induces a correlation, 

(5(^)S(k 2 ))\ hp{K) = / p (k 1 ,k 2 )/i;(K)e?.A*AJ^' iasI (1) 

where SP is shorthand for a Kronecker delta that sets 

K123 

ki + k 2 + K = 0. Note that h p (K) here are the new- 
field Fourier components during inflation when their ef- 
fect on primordial perturbations is imprinted. The func- 
tion / p (ki, k 2 ) is related to the density-density-new-field 
bispectrum B p (ki,k2,K) and new-field power spectrum 
P p (K) through B p (k u k 2 ,K) = P p (^)/ p (k 1 ,k 2 )e?.feifei 
Global SI requires that / p (k 1 ,k 2 ) be a function only of 
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FIG. 1: left six: The six possible polarizations to an otherwise statistically isotropic two-point correlation function. The first two 
modes are the usual transverse traceless tensor polarizations (gravitational waves), in which there are quadrupolar distortions in 
the plane transverse to the direction z of the wave. The next two are scalar and longitudinal- vector distortions, respectively. The 
scalar mode represents an isotropic modulation while the longitudinal-vector mode stretches and compresses the correlations 
along z. The two transverse-vector modes induce quadrupolar distortions in the xz and yz directions, respectively, right two: 
The circular polarizations of the tensor mode (h+t) and vector mode (h+ v ). 



kf, fcf, and ki • k 2 . 

The parameter p labels the polarization state of the 
new field and e^-(K) its polarization tensor, a symmetric 
3x3 tensor. The most general such tensor can be decom- 
posed into 6 orthogonal polarization states [H] , which we 
label s = {+, x,0,z,x,y}, that satisfy e^e p ^ = 28 pp i. 
These states can be taken to be two scalar modes e° ex Sij 
and ef • ex KiKj — K 2 /3, two vector modes e^- y ex K^w ' 
with K' l w l = 0, and two transverse traceless modes (the 
"tensor" modes) e + and e x . 

If K is taken to be in the z direction, then the + po- 



larization of the tensor mode has e+, = 
all other components zero, 



XX 
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polarization has 
£* = 1 with all other components zero. These two 
tensor modes are thus characterized by a cos 2<f> or sin 2<f> 
dependence, for e + and e x , respectively, on the azimuthal 
angle about the K direction of the tensor mode. The first 
two columns in Fig. [T] show the distortions induced to 
an otherwise isotropic two-point correlation function by 
correlation of the density field with a + and x polarized 
tensor mode. Shown there is a quadrupolar distortion in 
the x-y plane that then oscillates in phase as we move 
along the direction z of the Fourier mode. 

The first scalar mode has e?- = ^2/2>8ij and as shown 
in Fig. [IJ represents an isotropic modulation of the cor- 
relation function as we move along the direction z of the 
Fourier wavevector. The other scalar (or longitudinal- 
vector) mode has efj oc diag(— 1, — 1, 2)/-\/3 which repre- 
sents a stretching and compression along z. Both scalar 
modes represent local distortions of the two-point func- 



tion that have azimuthal symmetry about K. 

Finally, the two transverse-vector modes have e xz = 
e x zx = 1 with all other components zero, and e| z = e v zy = 
1 with all other components zero. These two modes rep- 
resent stretching in the ±xz and ±j/z directions, respec- 
tively, as shown in the last h x and h y columns in Fig. [T] 
These two transverse- vector modes have coscj) and sin <fi 
dependences on the azimuthal angle <f> about the direc- 
tion of the Fourier mode. 

The specific functional form of / p (ki,k2) depends on 
the coupling of the new field (scalar, vector, or ten- 
sor) to the inflaton. Global SI requires, though, that 
/ p (k!,k 2 ) will be the same for the two tensor polariza- 
tions and the same for the two vector polarizations; i.e., 
/ x (ki,k 2 ) = /+(ki,k 2 ), and / x (k!,k 2 ) = / y (ki,k 2 ). 
The same is not necessarily true for the scalar pertur- 
bations. In fact, the polar-angle dependence that distin- 
guishes the and z polarizations can be absorbed into 
/o(ki,k 2 ) and / 2 (k 1; k 2 ). Thus, in practice, one can de- 
scribe the most general scalar distortions to clustering in 
terms of either the or the z polarization by appropriate 
definition of /o(ki, k 2 ) or / z (ki, k 2 ). (This is the mixing 
between a scalar mode and a longitudinal- vector mode.) 
We thus below merge these two polarizations into a single 
polarization which we label with a subscript s. 

Suppose now that a correlation such as that in Eq. (nl) , 
for either a scalar, vector, or tensor distortion, is hypoth- 
esized. How would we go about measuring it? According 
to Eq. (IT]), each pair <5(ki) and <5(k 2 ) of density modes 
with K = ki + k 2 (note that we have re-defined the sign 



of K here) provides an estimator, 



mator, 



h p (K) = <J(k 1 )<J(k 2 ) / p (k 1 ,k 2 )ef,fcifc^ 



(2) 



for the Fourier-polarization amplitude h p (K). Since 
(|<5(k)| 2 ) = VP tot (k), where P tot (k) = P{k) + P n (k) is 
the measured matter power spectrum, including the sig- 
nal P(k) and noise P n (fc), the variance of this estimator 

is 



2VP tot (k 1 )P tot (k 2 ) / p (k lj k 2 )e?,.fci^ 



(3) 



The minimum-variance estimator for h p (K) is then ob- 
tained by summing over all these individual (ki, k 2 ) pairs 
with inverse-variance weighting: 



.— -~_ ^ /*(k, K - k)ef,-fc l (Zf - kV 

p{ ' p [ ' 2-. 2VP^(k)P^(\K - k|) 

x 6(k)d(K - k), 
where the noise power spectrum, 



Pp(K) 



2^ 2yp tot (fc)P tot (l K - k l) 



(4) 



(5) 



is the variance with which /i p (K) is measured. This 
P p (K) is a function only of the magnitude K (not its 
orientation) as a consequence of global SI, and for the 
same reason, P X (K) = P+{K) = P t (K), for both the 
signal and noise power spectra, and similarly P X (K) = 
P y (K)^P v (K). 

In general, the amplitudes /i p (K) arise as realizations 
of random fields with power spectra Ph{K) — AhP[ (K), 
for h — {s,v,t}, which we write in terms of amplitudes 
Ah and fiducial power spectra P h (K) . We now proceed 
to write the optimal estimator for the amplitudes A^. 

Each Fourier-mode estimator h p (K) for the appropri- 
ate polarizations (s for scalar, x and y for vector, and + 
and x for tensor) provides an estimator, 
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K,p 



Pl(K) 



y- 



h p (K) ~p;{k) 



(6) 



for the appropriate power-spectrum amplitude. Here we 
have subtracted out the noise contribution to unbias the 
estimator. If h p (J£) is estimated from a large number 
of (5(ki)-<5(k 2 ) pairs, then it is close to being a Gaus- 
sian variable. If so, then the variance of the estimator in 
Eq. |6| is, under the null hypothesis, 



P[(K) " [P;(K)\ 



(7) 



Adding the estimators from each Fourier mode with 
inverse-variance weighting leads us to the optimal esti- 



it P 2 [ p p( K )] 



v 



h p (K) -p;{k) 



where 



^ 2 =E[ p /w] 2 / 2 [w] 



K,p 



(8) 
(9) 



For the vector-power-spectrum amplitude A v we sum 
over p — {x, y} and for the tensor-power-spectrum am- 
plitude At over p = {+, x}. Following the discussion 
above, the sum on p is only for p = s for A s . 

The estimator in Eq. ([8]), along with the quadratic 
minimum- variance estimator in Eq. Q, demonstrates 
that the correlation of density perturbations with an un- 
seen scalar, vector, or tensor perturbation appears in the 
density field as a nontrivial four-point correlation func- 
tion, or trispectrum. The dependence of the trispectrum 
on the azimuthal angle about the diagonal of the Fourier- 
space quadrilateral distinguishes the shape dependences 
of the trispectra for scalar, vector, and tensor modes. To 
specify this trispectrum more precisely, though, requires 
inclusion of the additional contribution induced by modes 
K that involve the other two diagonals of the quadrilat- 
eral. Likewise, if a signal is detected — i.e., if the null- 
hypothesis estimators above are found to depart at > 3<r 
from the null hypothesis — then the optimal measurement 
and characterization of the trispectrum requires modifi- 
cation of the null- hypothesis estimators in a manner anal- 
ogous to weak-lensing estimators |13) . 

We now evaluate the smallest amplitudes A s , A v , and 
At that can be detected with a given survey. To do 
so, we take for our fiducial models nearly scale-invariant 
spectra P h (K) = A h K nh ~ 3 , with \n h \ < 1. We more- 
over take the density-density-new-field bispectrum to be 
of the form in Ref. [TJ. We then find that the inte- 
grand (using J^k ""*• V J d 3 k/(2n) 3 ) in Eq. ^ is dom- 
inated by the squeezed limit (K <C k\ — fc 2 ) where 
/ p (k 1; k 2 ) ~ — (3/2)P(fci)/fef. We then approximate 
P{k)/P tot {k) ~ 1 for k < /c max , where fc max is the largest 
wavenumber for which the power spectrum can be mea- 
sured with high signal to noise, and P(k)/P tot (k) ~ 
for k > fc max . This then yields a noise power spectrum 
P^ t} (K) ~ 207r 2 /£4 ax and P?{K) ~ 87r 2 /fcL x . Evalu- 
ating the integral in Eq. ([9]), we find the scalar, vector, 
and tensor amplitudes detectable at > 3c (for nu — 0) 
to be 



3(Th — 307rv 3irCh 



h'n 



288C h 



(10) 
where Cs t}V \ = 1 and C s = 2/5. The smallest detectable 
power-spectra amplitudes are thus inversely proportional 
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distribution can be decomposed into scalar, vector, and 
tensor distortions. We have presented straightforward 
recipes for measuring these distortions. Such effects may 
arise if the inflaton is coupled to some new field during 
inflation. We have avoided discussion of specific models, 
but the introduction of new fields during inflation is quite 
generic to inflationary models. We therefore advocate 
measurement of these correlations with galaxy surveys, 
and in the future with 21-cm surveys, as a simple and 
general probe of new inflationary physics. 

This work was supported by DoE DE-FG03-92- 
ER40701 and NASA NNX12AE86G. 



FIG. 2: The smallest scalar, vector, and tensor power- 
sepctrum amplitudes A s , A v and A t , respectively, detectable 
at the 3cr level as a function of the maximum wavenumber 
fcmax of the survey. Shown are results for survey volumes 
of 10 [Gpc/h] 3 and 200 [Gpc/h] 3 , or minimum wavenumbers 
fc min ~ 0.001 [ft/Mpc] and fc min ~ 0.003 [h/Mpc], respectively. 



to the number of Fourier modes in the survey. We show 
the projected detection sensitivities for surveys with vol- 
umes of 200 [Gpc/h] 3 and 10 [Gpc/h] 3 in Fig. § 

For example, if we apply this estimate to a tensor field 
and assume this tensor field to be primordial gravita- 
tional waves, then a sensitivity to a tensor amplitude 
A t ~ 2 x 10~ 9 near the current upper limit requires 
fcmax/fcmin ^ 5200. Such a dynamic range is probably 
beyond the reach of galaxy surveys, but it may be within 
reach of the 21-cm probes of neutral hydrogen during the 
dark ages envisioned in Refs. [TU1 [14] . Of course, the sig- 
nal could be larger if the inflaton is correlated with a 
scalar, vector, or tensor field that leaves no other trace. 

Finally, several new tests for parity-violating early- 
Universe physics can be developed from simple modifi- 
cation of the estimators above. To do so, we substitute 
the x and y polarizations, and + and x polarizations, 



with circular-polarization tensors e; 



±0 



e ij ± ie ij 



and 
e-j ± ief- . The two right-most patterns shown in 
Fig. [T] are the circular polarization patterns for tensor 
and vector modes. It may then be tested whether the 
power spectra for right- and left-circular polarizations are 
equal. For example, chiral-gravity models [15] may pre- 
dict such parity-violating signatures in primordial gravi- 
tational waves, and similar models with parity-violating 
vector perturbations are easily imaginable. 

Of course, "real-world" effects like redshift-space dis- 
tortions, biasing, and nonlinear evolution, must be taken 
into account before the estimators written above can be 
implemented, but there are well-developed techniques to 
deal with these issues [TB]- 

In summary, we have shown that the most general 
two-point correlation functions for the cosmological mass 
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